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We use an effective field theory (EFT) approach to calculate the next-to-leading order (NLO) 
gravitational spin-orbit interaction between two spinning compact objects. The NLO spin-orbit 
interaction provides the most computationally complex sector of the NLO spin effects, previously 
derived within the EFT approach. In particular, it requires the inclusion of non- stationary cubic 
self-gravitational interaction, as well as the implementation of a spin supplementary condition (SSC) 
at higher orders. The EFT calculation is carried out in terms of the nonrelativistic gravitational 
field parametrization, making the calculation more efficient with no need to rely on automated 
computations, and illustrating the coupling hierarchy of the different gravitational field components 
to the spin and mass sources. Finally, we show explicitly how to relate the EFT derived spin results 
to the canonical results obtained with the Arnowitt-Deser-Misner (ADM) Hamiltonian formalism. 
This is done using noncanonical transformations, required due to the implementation of covariant 
SSC, as well as canonical transformations at the level of the Hamiltonian, with no need to resort to 
the equations of motion or the Dirac brackets. 
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I. INTRODUCTION 



For some years now major efforts have been undertaken in order to detect gravitational radiation, one of Einstein's 
most important predictions of the theory of general relativity (GR), longstanding since 1916. Currently, three ground- 
based gravitational wave (GW) detectors are operating worldwide: LIGO (Laser Interferometer Gravitational-wave 
Observatory), Virgo, and GEO600 More advanced ground-based interferometers such as AIGO (Australian 

International Gravitational Observatory), that will complement the existing network of detectors due to its location 
in Australia or LCGT (Large-scale Cryogenic Gravitational- wave Telescope) in Japan Q, are expected to start 
operating in the near future. Moreover, LISA (Laser Interferometer Space Antenna) is expected to be the first 
space-based GW detector, exploring a complementary range of signal frequencies [g!]. 

One of the most promising sources of GW are binary systems of compact objects, in particular of black holes or 
neutron stars, that are found in the inspiral phase of their evolution. At this stage the dynamics of the binary can 
be described analytically in terms of the post- Newtonian (PN) approximation of GR ■ In order to succeed in the 
detection of GW signals, which is performed by matched filtering, accurate theoretical waveform templates of the 
inspiral are required, and so PN corrections should be obtained to high order, at least up to 3PN. The reason is that 
it is the overall value of the contribution of the PN correction to the phase of the waveform that matters for detection, 
and not its relative value, which can be numerically small with respect to the leading contribution. Since astrophysical 
objects are expected to have a spin, gravitational spin effects must be considered in the description of binary dynamics. 
Indeed, the leading order (LO) gravitational spin effect in compact binaries is evaluated to contribute already at 1.5PN 
[1, Q . Spin effects may enhance GW signals and increase the event rates for their detection [lH ■ 

A substantial possibility for progress in the analytical treatment was given when a novel effective field theory 
(EFT) approach for treating the PN formalism of GR was introduced by Goldberger and Rothstein [l^. The EFT 
approach is very advantageous in applying the efficient standard tools of quantum field theory to GR, notably handling 
the regularization required for higher order corrections in the PN approximation with the standard renormalization 
scheme. Moreover, this approach is appropriate to handle various problems that require the treatment of several 
typical length scales. Originally, the EFT approach was suggested for implementation on the inspiral phase of a 
binary to yield predictions of gravitational radiation [T2| . [l3| . Then, it was used to obtain the thermodynamics of 
higher dimensional Kaluza-Klein (KK) black holes and was extended for the rotating 15] and charged ^6i] cases. 

Various applications of the EFT approach then followed. Next-to-leading order (NLO) PN spinl-spin2 and spin- 
squared effects in inspiralling compact binaries were calculated for the first time ll7Hl9ll (the complete results were also 
obtained with the Arnowitt-Deser-Misner (ADM) Hamiltonian approach in l^-ljl). The 2PN dynamics for binaries 
were reproduced [23 |. and new results for n-body systems were obtained with automated calculations i_25]. The 
possibility to extract the three- and four-graviton vertices from binary pulsars and coalescing binaries was examined 
in [2^. Radiative corrections in compact binaries were reproduced up to 3PN order 27 1. Self- force on extreme mass 
ratio inspirals and radiation reaction were also treated in the EFT approach \2§, fSoj . Most recently even finite size 
corrections to the radiation reaction in classical electrodynamics were found using the EFT approach 30]. 

In this work, we make an application of the EFT approach to the calculation of the NLO spin-orbit (SO) interaction 
in inspiralling compact binaries. This is the only NLO PN spin correction that was not handled before in the EFT 
approach due to its higher complexity. In particular, the NLO SO sector requires the inclusion of nonstationary 
cubic self-gravitational interaction, as well as the implementation of a spin supplementary condition (SSC) at higher 
orders. Here, we work with the reduction over the time dimension that was proposed within the EFT approach by 
[3H , and use the nonrelativistic parametrization of the gravitational field, which was first utilized in the computation 
of the NLO spinl-spin2 interaction [l^. The nonrelativistic gravitational (NRG) field variables prove to greatly 
simplify the calculation and make it feasible with no need for automated computations. The NRG fields also illustrate 
the coupling hierarchy of the different gravitational field components to the spin and mass sources. We make an 
effort and put an emphasis on presenting the treatment of self-gravitating spinning objects in an EFT approach 
in a clear and self-contained manner. Then, after we give a detailed account of the EFT calculation of the NLO 
SO sector, we use noncanonical transformations, required due to the implementation of covariant SSC, as well as 
canonical transformations at the level of the Hamiltonian. This leads us to the canonical result calculated in the 
ADM Hamiltonian formalism ^] , which also agrees with that of [33] : whose authors first obtained the complete NLO 
SO interaction equations of motion (EOM). Naturally, it is more efficient to use variable transformations and remain 
in the level of the Lagrangian/Routhian, or Hamiltonian, than to resort to the treatment of EOM and the Dirac 
brackets. 

Throughout this paper, we use c = 1, t]^^ = diag[l, —I, ~l, ^l], and the convention for the Riemann tensor 
is R'^^ajj = da^'^jj — dpV^^ -\- r^^r^^ — r^^r^^. Greek letters will denote indices in the global coordinate frame. 
Lowercase Latin letters from the beginning of the alphabet will denote indices in the local Lorentz frame, while 
uppercase Latin letters denote the body fixed, corotating Lorentz frame. All these indices will run over to 3, while 
spatial tensor indices over 1 to 3, will be denoted with lowercase Latin letters from the middle of the alphabet. The 
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letter t may be used as an alternate for the time index 0. We will also use the following notation = J ^^J^^ for 

abbreviation (boldface characters are used to denote 3- vectors). 

The paper is organized as follows. In Sec.|ll]we review the EFT approach for the binary inspiral problem in the PN 
approximation, with the NRG field parametrization, and present the Feynman rules in the linearized harmonic gauge, 
required for the EFT computation. In Sec. Illll we present in a self-contained manner the formalism for the treatment 
of self-gravitating spinning objects in an EFT approach and present the Feynman rules involving the spin required 
for the computation with the NRG fields. In Sec. IIVI we derive the LO SO potential, which is required to obtain the 
NLO SO Hamiltonian, but also to illustrate the most simple EFT calculation in the SO sector and its ambiguities. 
In Sec. |V]we present in detail the full calculation of the NLO SO Lagrangian/Routhian, going over all contributing 
Feynman diagrams, from one- and two-graviton exchanges to cubic self-gravitational interaction diagrams. In Sec. IVII 
we present the NLO SO Lagrangian/Routhian EFT result and derive the NLO SO Hamiltonian, using noncanonical as 
well as canonical transformations to obtain the canonical result in [32] in the ADM gauge. In Sec. IVIII we summarize 
our main conclusions. 

We note that after this work was complete, there appeared [s^, HI] , which study the same problem. 



II. EFT APPROACH FOR BINARY INSPIRAL IN THE PN APPROXIMATION 

As explained in the EFT description of the binary inspiral problem [l^ . expanding the metric around flat spacetime, 
the metric can be decomposed into potential and radiation modes, both with the same typical time variation scale of 
the binary orbital frequency v/r, but with a different typical length scale for each, r and r/v, respectively, where v 
and r are the typical orbital parameters of the binary, and we are working in the limit w <C 1 of the PN approximation. 
Thus, as far as the potential gravitons are concerned, the radiation component of the field is just a slowly varying 
background field of soft momenta gravitons. Moreover, considering these typical scales, we see that the potential 
gravitons are off shell with their frequency being much smaller than their momentum, and hence can be approximated 
as stationary to leading order. This motivates a KK reduction over the time dimension for the potential field modes 
as suggested in [3, [3l| . Note that since we are concerned with conservative dynamics in this work, we only need the 
potential modes anyway and can simply set the radiation modes to zero. 

Therefore, instead of the common Lorentz covariant metric parametrization in the weak field limit, of the form 
gfiii — Vfii' + ^M'' ' metric is parametrized in a nonrelativistic form according to the Kaluza-Klein ansatz 

dr^ = g^^dx^^dx" EE e2^(dt - A, dx'f - e-'^'''-f^jdx'dx^ . (1) 

This defines a set of new fields (0, Ai,"/ij), with a scalar field corresponding to the Newtonian potential, a 3-vector 
field Ai, corresponding to the so-called gravito-magnetic vector, and a 3-dimensional symmetric 2-tensor field jij, the 
nonrelativistic gravitational fields as discussed in [3l[ . Then, in terms of the NRG fields the metric reads 

\ -e^^A, -e-^H, + e'^A.A, )-\-A,- 2A,(j) -S,, + 2^S,, - a,, - 2cf>^S,, + 2^a,, + A,A, J ' 

where jij = Sij + Uij , and we have written the approximation for the metric in the weak field limit up to second order 
in the fields, which is the order required in this work. In the ground state (0) = {Ai) = (atj) = 0, and the metric 
reduces to Minkowski spacetime. We should also have the form of the inverse metric, and it is given by 



^ , - e^^-f'^AA, -e^^A^ \ / 1 - 20 + 20^ _ AkAk -A, - 2Mj + ^jkAk 

^ \ -e^'^A' -e^'l'fi ) -\ -A, - 20yl, + a.kAk - 20% -f fT„- - 202^^- + 20(t„- - aaCTi] 

(3) 

where Y^jjk = Sik, i.e. 7*^ is the inverse of jij, and A' = j''^ Aj. Here too, we have written the weak field limit of the 
inverse metric up to second order in the fields. 

Next, we consider the purely gravitational action. It is the usual Einstein- Hilbert (EH) action plus a gauge fixing 
term of our choice 

Sg Seh + Sgf = - ^Q^Q J '^'^^ iVsR + ^gf) ■ (4) 

We suppress the time dependence of the fields in the EH action to obtain the KK reduced action for the gravitational 
field, given by 



SkK = SEH\stat. — 

I67 



i- / dtd^x V7 \-R[l^,] + 27''' 5.^-0 - ]:e^^F,,Fuf''l'' 



(5) 
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where 7 = det{'jij) and Fij = diAj — djAi. 

Now, we must include a gauge fixing term to set the gravitational action completely. We should be able to invert 
the graviton kinetic term, i.e. the quadratic part of the Lagrangian, in order to obtain the graviton propagators. 
Then, the propagators can be read off from the terms that are quadratic in the fields, depending on the choice of 
gauge fixing term. Quadratic terms with time derivatives are suppressed as subleading corrections in powers of u^: 
We recall that the potential gravitons are instantaneous within the leading stationary approximation, representing 
off shell gravitons. A time derivative adds a factor of v for potential modes, hence the propagator corrections are 
suppressed by v'^. The quadratic part of the KK reduced action (after two integrations by parts) is given by 



S 



1 



KK[g^ 



327rG' 



dtd X 



(6) 



Hence, we take the natural gauge adequate for the NRG fields, namely, the Lorentz gauge for the vector field Ai, and 
the harmonic gauge for the 3-dimensional 2-tensor field aij. This is just equivalent to the harmonic gauge for the 
Lorentz covariant parametrization of the metric. Thus, the gauge fixing term is given by 



5, 



1 



GF[KK] 



32nG 



dtd^x 



{^^A,Y 



(7) 



Therefore, the NRG scalar, vector, and 2-tcnsor field propagators in the harmonic gauge are given by 

3ik-(xi-X2) 



(a;i) (I3{x2) ) 



{A,{xi) A,(x2)) = -- 



1 
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1 
2 



= {aij{xi)aki{x2)) = Pir,k 



„lk-(xi-X2) 

^ik-(xi-X2) 
k2 



5{ti-t2), 
S{tl-t2), 
S{tl~t2), 



(8) 
(9) 
(10) 



where Pij-ki = \ (SikSji + Sudjk — 2SijSki)- Here and henceforth, the Feynman rules are presented in position space. 
Note the simple form of the propagators obtained with the NRG field parametrization and harmonic gauge, in 
particular, that of the (j) and Ai fields, which dominate in the interaction, an advantageous feature of the NRG field 
variables. 

However, as already noted there are corrections to the nonrelativistic instantaneous graviton propagators, i.e. to the 
graviton kinetic terms in the KK reduced action. These are included systematically as perturbative corrections, where 
each propagator correction is suppressed by a power of v for each time derivative it contains. It is a departure from 
the stationarity approximation in the gravitational action, thus where the use of the KK action becomes insufficient, 
and we are bound to go back and consider the full EH action. Extracting from the EH action the quadratic time 
dependent terms, which contribute to the order considered here, namely, the terms which contain the scalar field (p 
as well as the vector field A, , we have 



1 



167rG 



6{do^f + 4d^A,do<l> + doA,d,<j^j 



^doAidiajj - ^diAidoCTjj - 2do(j)docri., 



(11) 



Thus in order to simplify as much as possible, we refine the gauge fixing term of Eq. ([7]) to eliminate undesired 
quadratic vertices. Hence, we reset it to 



1 



'GF 



327rG 



d^x 



d,A, + ( 4300 - \d^o. 



d^A 



(12) 



so that to the order considered here, we are left with the following time dependent quadratic terms for the scalar and 
vector fields 



1 



327rG 



d^x {(d^A.f - 4(ao</.)2 



(13) 



Again, this is just equivalent to the harmonic gauge for the Lorentz covariant parametrization of the metric. Note that 
here we face the major advantage of the NRG parametrization in harmonic gauge - the 2-point functions between the 
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three different fields are 0: {(j>Ai 
vertices are given by 



ajk) = {Ai<Tjk) = 0. Thus, the Feynman rules for the propagator correction 



1 



SttG 
1 



(14) 
(15) 



The encircled crosses represent the self-gravitational vertices that contain two time derivatives. 

We will also have contributions from three-graviton vertices of cubic gravitational self-interaction. These are 
extracted from the purely gravitational action as we are using the linearized harmonic gauge, and not the fully 
harmonic gauge, i.e. since our gauge fixing is quadratic in the fields, we extract the cubic part from the KK or EH 
action alone. The LO three-graviton vertices scale as according to the EFT power counting |12|], and they are 
easily derived from the cubic part of the KK reduced action. To the order considered here, there are two such cubic 
vertices contributing, the (j)F'^ and the cr0^ vertices given by 



1 



167rG 



dtd^x [(-20 {diAjd^Aj - diA-^djAi)) + {~2atjd^(j)dj(j) + a^jd^^ 



(16) 



However, at the order considered, we also have a contribution from another three-graviton vertex, including yet again 
a departure from the KK stationary action. This vertex contains a single time derivative, hence it scales as v^. 
Therefore, to obtain this vertex we have to go through the tedious extraction of the cubic part of the EH action, 
containing a very large number of terms, each with a complicated tensor index structure as was done in [l^ . The 
contributing time dependent cubic vertex is the Acj)^ vertex given by 



S 



EH 



1 



IGttG 



dtd^x [4Aid,(t)do(l)] 



Hence, the Feynman rules for the three-graviton vertices are given by 




1 



SttG 
1 

167rG 
1 



x)diAj{x) diAj{x) - djAi{x) 



d^x [2ai-i{x)di(j>{x)d-i(j>{x) - crjj{x)di(j){x)di<p{x)], 



d'^x [Ai{x)di(j)[x)d(i(j){x)]. 



(17) 

(18) 
(19) 
(20) 



The gray encircled cross represents the self-gravitational vertex, containing a single time derivative. 

Now we consider the gravitational coupling to the two massive compact objects. We take the worldline action of a 
point particle (PP) for each of the two objects, so that to the order considered here we have for the binary system 



pp{m) 



2 

E 

n=l 



ninddn 



(21) 



where the ellipsis denotes subleading operators encoding finite size effects. Here, we do not consider them since they 
enter only at 5PN |12|. We parametrize the worldline using the coordinate time t — x'^ ^ i.e. ct = i, so that we have for 
= dx^ /da: vP = 1, = dx'^/dt = w*. Expanding the PP action of a massive object using Eq. ([T]), up to the order 
required here in the fields and velocities, gives rise to the following couplings to the worldline mass: 



Spp{m)= - J mda = -m 
1 



dt 



dx n dx u 



dt 



dt dt 
3 



dt 



Ay + |</)t;2 



1 



(22) 



where the ellipsis stands for one-, two-, and n-graviton couplings, and velocity powers beyond the order considered 
here. It is clear that any number of gravitons can couple to the worldline mass. Moreover, the LO couplings of 0, 
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Ai, and (Jij are 0(w°), 0{v^), and O(w^), respectively, and all couplings contain an infinite power series in v"^. The 
Feynman rules for the one-graviton couplings to the worldline mass are given by 



m dt (^(x(t)) 



1 + -vit)' + 



m dt A,{-K{t))v\t) 



= m I dt a^■j{x{t))v'■{t)v^t) 



-Mtf 



(23) 



(24) 



(25) 



where the heavy solid lines represent the worldlines, and the spherical black blobs represent the masses on the 
worldline. Again, the ellipsis denotes higher orders in w, beyond the order considered here. For the two-graviton 
couplings to the worldline mass, we have the following Feynman rules: 



= -m / dt (t>{^{t)f 



m / dt (l)iTc{t))A^{ji{t)y(t) 



(26) 



(27) 



III. SELF-GRAVITATING SPINNING OBJECTS WITH AN EFT APPROACH 

The interaction of classical spin with gravity was discovered in 1937 by Mathisson [s^. However, the EOM for a 



spinning test body in GR were derived in their modern form by Papapetrou [37[. The Mathisson-Papapetrou (MP) 
EOM were obtained based on a multipole formalism for extended bodies and describe the motion of a spinning body 
in the so-called pole-dipole approximation, where quadrupole and higher multipole moments are neglected. Later, the 
MP equations were rederived using an action approach. For example. Bailey and Israel obtained the MP equations 
from a variation of the action, though they did not speciiy the Lagrangian form [s^ , similar to the analogous treatment 
given in flat spacetime by Hanson and Regge Yee and Bander used a Routhian, whose form they specified, to 
obtain the EOM fio^ . More recently, further similar treatments using an action approach were made in |l3, 14T1 - |43| . 
Here, we consider the self-gravitating spinning objects of the binary system as point particles in what is equivalent to 
the pole-dipole approximation, i.e. we only work to linear order in the particles spins. That is ,we are not concerned 
with finite size effects, which are quadratic in the spin, such as self-induced effects or tidal deformations, or in any 
intrinsic permanent multipole moment beyond spin. 

We start by defining a tetrad field as a set consisting of a timelike future-oriented vector and three spacelike vectors, 
which satisfy 



5M-e^^e^ = r,ab, (28) 
with rjab = diag[l, —I, —1, —1] the fiat space Minkowski metric, and the set of vectors with the tetrad index a 



ranging from to 3. Tetrad indices are raised and lowered with the Minkowski metric rjab, e.g. e"^^ = ry"''e^'. This is 



equivalent to defining the reciprocal tetrad field by 

g^.e'^^el = (29) 
which also tells us that the tetrad is an orthonormal set. Now, one can easily obtain the completeness relation 

Vabe'^'el = Si;. (30) 

From Eqs. (I28l) - (p0)) it is easy to obtain the following expressions for the metric and its inverse in terms of the tetrad 
fields: 

9^. = Vabelel g^'^ = rf''e^^el. (31) 
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The projections of a vector (and similarly for a tensor of any rank) onto the internal tetrad space are defined as 
ya ^ gayM^ ^ eJ^V^, or convcrscly we have Vf" = e^F", V,, = e^Va- 

In order to describe the rotation of a spinning object in curved spacetime, we need two different tetrads. First, we 
need the tetrad of the body fixed, corotating Lorentz frame, carried by the rotating particle. To begin with we would 
like to follow the evolution of this tetrad along the worldline. We would also require a reference tetrad that covers 
the background spacetime, the tetrad of the local Lorentz frame of the particle. The internal tetrad space is Lorentz 
invariant, so while satisfying Eq. (PT|) . we have a freedom in the choice of this tetrad that will be fixed. We recall 
that our convention for the notation of the body fixed tetrad is by Latin uppercase letters from the beginning of the 
alphabet, whereas the background tetrad is denoted by Latin lowercase letters from the beginning of the alphabet. 
The body fixed tetrad is related to the background tetrad by a Lorentz transformation, i.e. = A^e^, where the 
Lorentz transformations and their inverse are defined by K^°'K^^riAB — rj'^'' and KAa^BbTj'^^ = rjAB- The Lorentz 
transformations relating the two tetrads should reduce to pure rotations in the local rest frame. Now we can proceed 
to define the generalized angular velocity in curved spacetime, which describes how the body fixed tetrad rotates 
along the worldline, by 

n^'^ ^ 4-^, (32) 

where D/ Da is the covariant derivative with respect to the worldline parameter a. It is easy to verify that fi^'' is an 
antisymmetric tensor. This definition is a generalization of the definition for flat spacetime given by VL"^^ = '^^^ 
[39]. 

We should construct an action for the spinning particle that respects general coordinate invariance, reparametriza- 
tion invariance (of the worldline parameter of the objects), as well as SO (3) residual Lorentz invariance of the 
tetrad internal space. The minimal coupling part of the action would then be given by the PP action Spp = 
J da Lpp(si, S2, S3, Si) = J da Lppiu'' T^^'^), where the Lagrangian Lpp is a function of the four covariant scalars 
Si, . . . ,S4 that can be formed from the coordinate velocity = dx^ /da and the generalized angular velocity il'^'^ 
given by [s^ 

S2 = nt^'^n^u, 

Si = nt'^n^piip'^n^p, or det{n'"'). (33) 

Since the action is reparametrization invariant, we require the Lagrangian Lpp to be a homogeneous function of degree 
1 in the velocities u'^ and fi'"^ [39*1 . So that if we define the linear momentum p^, and the antisymmetric tensor 5'^^ 
of the particle as 



we have from Euler's theorem that Lpp is given by 



dL 

duf' 
, dL 



(34) 
(35) 



Lpp = -p^m'' - -S^,,n^'r (36) 



Note that in curved spacetime the linear momentum is different than the momentum conjugate to the coordinate 
velocity m^, since $7'^" depends on (e.g. is explicit in the covariant differentiation in the definition of VL^^'^ in 
Eq. ([211)). Therefore, the differentiation in the definitions of Eqs. ([M]) and ([55]) is special in that we treat Vl^^'^ as 
independent of w^. The minus signs in these definitions are chosen to give the correct nonrelativistic limit. Indeed, 
the PP Lagrangian form of Eq. ([55)1 is just obtained from a minimal coupling of the flat spacetime action |39| . 

A variation of the action with respect to the position a;^, and the tetrad field e^, with fixed endpoints yields the 
EOM for the position and spin respectively [s^, |4l| . The EOM are given by the well-known MP equations 

^ = -Ir.up.W^S^^, (37) 



Da 



P^iUy-pyUp. (38) 
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However, the description of the motion of a spinning object by the MP equations is not closed. The sohition for 
the EOM is not unique due to redundant unphysical degrees of freedom associated with the spin that should be 
eliminated. Hence, we need to choose a so-called spin supplementary condition to close the description. The choice 
of SSC can be understood as the freedom to choose the point internal to the spinning body, whose motion along the 
worldline is followed, as any spinning "particle" must actually be an extended object, and has a lower bound for its 
finite size due to its spin. However, the SSC should fix the representative worldline of the body uniquely^ whereas not 
all choices of SSC in curved spacetime do [i^. Tulczyjew was the first to give a definition of SSC, i.e. of a center of 
mass, which ensures the uniqueness of the EOM [45]. 

An appropriate SSC in curved spacetime may be written in the form 

C = S'^'V^ = 0, (39) 

where is some suitably chosen timelike vector For example, Tulczyjew suggested the covariant SSC given by 

El 

S^-'pu = 0, (40) 



where p^, is the linear momentum here also. Note that only three out of the four components of Eq. ((39)) are 
independent constraints as C^V^ = 0. This corresponds exactly to the three redundant spin degrees of freedom that 
should be eliminated. Thus, the three independent constraints are taken to be the spatial components C\ Also, is 
a zero vector, so of course the SSC can be applied in any preferred frame. In particular we may choose to apply it in 
the local Lorentz frame, or in the coordinate frame. The preservation of the SSC upon evolution, i.e. the requirement 
that DC^/Da = 0, together with the MP equations Eqs. dST]) and applied to the covariant SSC in Eq. (j40|), 

yields the linear momentum p'^ in terms of the coordinate velocity and the spin 

p'^ =m—+0{S^), (41) 



where the part quadratic in the particle's spin, not written explicitly, is SSC dependent. We see that if we plug p'^ 
from Eq. (PT|) into Eq. (1551) . we recover the nonspinning part of the PP action given in Eq. ([2^ . Also if we plug 
Eq. (j4T|) into the covariant SSC in Eq. (|40|) . we see that in our approximation the covariant SSC can be taken as 

S^"'u^ = 0. (42) 

Any SSC of the form of Eq. ([39)1 . which guarantees the uniqueness of the solution of the EOM, can be used 
to eliminate the redundant degrees of freedom of spin, namely, in the end to be substituted in the temporal spin 
components Sot. The basic procedure for this to be carried out is to add a Lagrange multiplier vector multiplied 
by the SSC, i.e. a term A^C, to the Lagrangian, as was done in e.g. However, the use of Lagrange multipliers 
to incorporate the SSC in the action would be out of place in a treatment of the problem in an EFT approach. Yee 
and Bander, who formally worked with a Routhian due to the explicit use of the spin momentum as an independent 
variable, with its conjugate coordinate being cyclic [i^, added an acceleration dependent piece to the Routhian given 
by 

U^-Sab Dub , - 

in order to ensure that the covariant SSC of Eq. (|42|) is conserved fioj. Note that this term is proportional to the 
SSC itself. Following the work of Yee and Bander a new algebraic Riemann dependent operator was also added to a 
Routhian in p7| . The additional operator given by 



is obtained from the term in Eq. (|431) . by using the MP equation in Eq. (|37| . and the relation in Eq. (j41l) . This term 
is quadratic in the spin of the particle, i.e. of the form of a self-induced quadrupolelike operator. However, the main 
disadvantage of working similarly with the SSC in Eq. (|42l) is that it leads to noncanonical variables, i.e. such that the 
EOM, in particular, the spin EOM, which are conveniently obtained from a Hamiltonian/Routhian, cannot be obtained 
by using the usual canonical Poisson brackets, and one will have to make some noncanonical variable transformations 
to arrive to canonical variables that satisfy the canonical algebra. Alternatively, the SSC can be imposed within the 
Hamiltonian formalism by replacing the original Poisson brackets with the so-called Dirac brackets [39] . In that case 
the constraints in Eq. (j39p alone are not sufficient to describe the phase space consistently, as one needs to impose an 
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additional set of three constraints on the configuration space variables conjugate to the spin, so that the constraint 
hypersurface in phase space contains the same number of configuration coordinates and conjugate momenta. The 
Dirac brackets are in essence the projection of the original symplectic structure onto the phase space hypersurface 
defined by the constraints. We note that in an action approach, there is no need for additional conjugate constraints, 
as was made clear in [42l] for example, where the corresponding Lagrange multipliers are found to vanish, i.e. these 
constraints are redundant. Finally, the constrained Hamiltonian will be obtained simply by inserting the constraints 
directly in the original Hamiltonian. 

The so-called Newton- Wigner (NW) SSC is well known to be the only SSC in flat spacetime, which yields canonical 
variables [H, HI] . Recently, the NW SSC was generalized from flat to curved spacetime in [H, |4§| . In [i^ they 
show that their generalized NW SSC leads to canonical variables at linear order in the particle's spin. A generalized 
NW SSC may be formulated as 

S^-'{p,+me°)^0, (45) 

where the tetrad here is the local background tetrad, and Pi, is again the linear momentum. Note that Eq. (P5|) is of 
the form given in Eq. (p9|. and that clearly Eq. (145]) reduces to the fiat spacetime NW SSC given by 



S'^^'p^ + mS''° = 0. (46) 

It should be noted however that a generalization of the NW SSC to curved spacetime is not unique. Here though, 
we will work with the covariant SSC of Eq. (|42l) . which are validly incorporated into our action approach. In the 
covariant SSC of Eq. , we have for the spin temporal components that 

S^O ^ s^3^J^ (47) 

which implies a suppression of the temporal spin components S''^ by one order of v with respect to the spatial 
components S^^ . We note that here we apply the SSC in the local Lorentz frame, since we will work with the spin 
variable that is projected onto the local tetrad space, so the velocity here is the local velocity, and not the coordinate 
velocity. Similarly, to the order considered here in the generalized NW SSC of Eq. ((45|) we would have 



S^o ^ K^J^J ^ Is'^v^v'^ + ■■■ , (48) 
2 8 

which implies a suppression of the temporal spin components S"* by at least one order of v with respect to the spatial 
components S^^ , and there is an infinite expansion in the local velocity v. 

Let us consider then the power counting of spin within the EFT description. The power counting rules in the EFT 
description of a binary that does not have spinning components within the PN approximation were given in (l^ IT3| , 
and we will not repeat them here. However, we should address here the power counting of a spinning compact object, 
in which we are concerned. We note that the spin of a compact object scales as 

S ~ mVrotTs < mrs ~ mv^r ~ Lu, (49) 

where Vrot is the rotational velocity of the object, Ts ~ Gm is its Schwarzschild radius, and we have used the Virial 
theorem Gm/r ~ . We therefore see that the spin gets suppressed by one order of v with respect to the orbital 
angular momentum, even for maximally rotating objects {vrot — 1), which we will take our objects to be. However, 
there exists another way of PN counting the spin, the so-called "formal" counting |49l] . In the formal counting the 
PN orders are counted from the original Einstein equations. Our counting adds half a PN order for each spin variable 
appearing in the expressions, compared to the formal counting. For example, the NLO SO and spinl-spin2 sectors are 
of the orders 2.5PN and 3PN, respectively, in our counting, but are both of the 2PN order in the formal counting. Our 
counting refiects better the true relevance of the PN spin corrections, and does not overestimate their importance. 
Moreover, it corresponds to the standard practice when defining the templates of LIGO/ Virgo and LISA detectors. 
However, it does not refiect well the computational complexity involved. In the formal counting some similarities to 
PN calculations for nonspinning objects are more manifest. To some extent, the formal counting refiects better the 
computational demands, e.g. the difficulty of the integrations required, and the regularization techniques that must 
be employed. 

Finally, let us focus our attention back on the linear in spin part of the PP Lagrangian given by the second term 
in Eq. (|36p . It can be rewritten in terms of the Ricci rotation coefficients (also called the spin connection) defined by 
(up to a sign convention) [Hoj 



(50) 
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The relevant spin part of the PP Lagrangian ipp(s) then reads 

LppiS) = -ISati^''' - i^ab<^^ (51) 

We recognize that the first term here represents the rotations of the body fixed Lorentz frame with respect to the 
background Lorentz frame. Hence, this kinetic term of flat spacetime, which docs not contribute to the interaction 
of spin with gravity, wih be omitted henceforth. Note that the spin that appears here is the projection of the spin 
tensor onto the background local tetrad, i.e. S'"'' = S^'^e'^e'^,. Therefore, the gravitational coupling to the spins of the 
objects is described by the relevant spin part of the PP action for each of the two objects in the binary system, so 
that we have for each [i^l 

SppiS) d^S'^'^f^abU^' + ■■■ . (52) 

where the ellipsis denotes subleading corrections beyond linear in spin, or due to finite size effects. Now, since it is 
convenient to obtain the EOM for the position of the particle from a Lagrangian, whereas the EOM for the spin are 
easily obtained from a Hamiltonian, we may consider the spin momentum as an independent variable, where we note 
that its conjugate coordinate is cyclic. Thus, we can consider our Lagrangian as a Routhian for the purposes of EOM 
derivation, and the Legendre transformation to a Hamiltonian with respect to the spin conjugate coordinate is trivial. 

In order to derive the Feynman rules for the spin couplings, it is convenient to use the common Lorentz covariant 
parametrization for the expansion of the metric around flat spacetime g^i/ = + /i^^. Then, using Eq. pip, one 
can flnd the general solution for the background tetrad, expanding in terms of h^^^ 



e;=K + ^ VA"^" - Ih^ph'.A^'' + ■■■■ (53) 



However, we can use the internal Lorentz invariance of the background tetrad to fix our reference tetrad to be aligned 
with the Lorentz frame at infinity so that 



el^S; + -hi - -h^hP, + ■■.. (54) 



2 

Using this tetrad, and expanding in Eq. (|52l) up to second order in h^^, which is all we will need here, we obtain the 
following Lagrangian: 

Going back to our NRG field parametrization using Eq. ([2]), and our choice for the worldline parameter cr = i, we 
obtain the following couplings to the worldline spin: 

+S'^Fij<j> - ^S'^A,dj(l) + 2S^'dM + • • • , (56) 

where the ellipsis stands for one-, two-, and n-graviton couplings, and velocity powers beyond the order considered 
here. Similarly to the mass couplings, any number of gravitons can couple to the worldline spin. Here however, unlike 
the couplings to the worldline mass, the order in the velocity v of the couplings is not explicit, due to the 5*°* spin 
entries and the freedom in the choice of SSC, in addition to the possibility of having time derivatives on the fields, 
since the spin is derivative-coupled. Also, here it is the Ai field that couples at LO in w, whereas the 4> and <Jij fields 
are subleading couplings, both at the same order of v (hence, the assignment "gravito-magnetic" for Ai, due to its 
dominant coupling to the spin). 
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The Feynman rules for the one-graviton couplings to the worldline spin are thus given by 



dt 



(57) 



= / dt 



(58) 



= / dt 



(59) 



where the gray oval blobs represent the spins on the worldlines. Note the time derivative appearing in the coupling of 
the Ai field to the worldline spin in Eq. (j57p . This is a departure from the leading stationary approximation. Note also 
the 5*°* spin entry that appears at the leading order coupling in Eq. ((55)) . This coupling enters already at the LO 
spin-orbit interaction (as we will see shortly in Sec. lIVp . so the use of SSC is required already at the LO gravitational 
spin correction at 1.5PN. Moreover, this S'°*0 coupling adds much complication to our calculation. We also note that 
the ellipsis in Eq. (|58p stands for another spin coupling term with a time derivative, which equals S^'^docfiv^. However, 
it can be seen from power counting that to the order considered here this contribution would eventually drop at the 
implementation of the SSC. 

For two-graviton couplings to the worldline spin the Feynman rules are 




= / dt 



= / dt 



2S^\t)d.,A,{^{tm^{t)) -5'M<)A,(x(t))a,</)(x(i)) 



2^•'X^)9^</'(x(^))0(x(^)) 



(60) 



(61) 



Note that the two-graviton ip(p coupling appears with the 5°' spin entry, i.e. the use of SSC is required here. Also 
recall that the velocity, appearing in the SSC derived expressions in Eqs. (|47|) and (|48|) for the 5*" entries, is the 
local velocity and not the coordinate velocity, appearing here in the Feynman rules. So we actually have metric field 
corrections in the temporal entries from the SSC rewritten as 



(62) 



For example, for the covariant SSC in Eq. (|47]), we have in terms of the coordinate velocity S'**' — S^^v^ - 
with the gravito-magnetic vector field Aj as the lea ding order metric field correction in the temporal spin entries, 
as was already found in the NLO spinl-spin2 sector 'l8||. These metric dependent corrections may be considered as 
contributions to the two-graviton couplings in Eqs. ()60p and (|61|) here, arising from the one-graviton spin coupling in 
Eq. (|58|) with 5*'° entries, if the SSC is applied at the level of the action. 

It is important to note that the addition of the algebraic Riemann dependent operator from [17], appearing in 
Eq. (j44|) . will not affect the Feynman rules relevant for the SO sector (nor for the spinl-spin2 sector), since the 
operator is quadratic in the spin of the particle. Finally, it should be pointed out that the fact that the spin is 
derivative-coupled also raises the complexity of the computations of spin corrections. 



IV. LEADING ORDER SPIN-ORBIT POTENTIAL 

In this section we derive the well-known LO spin-orbit potential, originally computed in [1, |^ . The result here will 
be important later in Sec. lVIl in the derivation of the NLO spin-orbit Hamiltonian. Moreover, we use this most simple 
computation of the SO sector to illustrate the computation of Feynman diagrams that contribute to this sector, and 
its ambiguities related with different choices of SSC. We will show how these ambiguities can be resolved in terms 
of canonical and noncanonical transformations. The SO sector contains the diagrams, which have one spin coupling 



12 



of one of the particles, together with at least one mass coupling of the other particle in the binary. First, we should 
use the power counting rules of [l^, [3], and our power counting for spin given in Eq. P^)) . to assign a power in the 
velocity v in the PN expansion for each term in the Feynman rules given in Secs.lllland lllll Then, we construct the SO 
diagrams, which contribute to a certain PN order, by considering all diagrams composed of all possible components, 
whose PN order adds up to the relevant PN order. The LO spin-orbit sector contains contributions only at 0{G), at 
which we have the single topology of one-graviton exchange. 

From our PN expansion of the Feynman rules, we find that the LO spin coupling is given by the first term of 
Eq. (HH), i.e. by 

i(s)2 = \s''F,j, (63) 

where the numerical subscript denotes the power in v of this term in the PN expanded Feynman rules. This can only 
be contracted, using Eq. with the LO mass coupling of Ai that is given by 

L(„,)i = mA,v\ (64) 

We obtain a one-graviton exchange diagram depicted in Fig. 1(a). It has the value 

Fig.l{a)^ _lSi • V2 X n+ [1 2], (65) 

where we denote r = xi(t) — X2(t), r = |r|, and n = ^. The spin is represented by a 3-vector defined by S"^^ = e^^^S^ . 
The labels 1 and 2 are used for the left and right worldlines, respectively. The [1 O 2] notation stands for a similar 
term, contributed by the mirror image of the diagram displayed, whose value is obtained under the interchange of 
particles labels. Note that under this exchange n — > — n. Although r, V2, Si, etc. depend on t, we suppress this 
dependence here and henceforth. Moreover, here and henceforth a multiplicative factor of J dt is suppressed and 
omitted from all diagram values. Here, we had to evaluate the Fourier integral coming from the propagator, and for 
that we have used Eq. (IA2|) in Appendix 1X1 It is easy to see that this diagram scales as (1.5PN). Hence, considering 
the next order in the spin couplings, we realize that there is another diagram contributing to the LO SO potential. 
We go on to the spin couplings of order given by Eqs. ((55)) and (ICTl) . so that we have 

L(s)3 - S'^d.^v' + + ]^S'^d,n,kv\ (66) 

where here the temporal spin entry S"* must also be taken into account, considering the SSC. However, contracting 
the aij field with the LO mass coupling of Uij given by 

L{m)2 = ^ma.jV^^ (67) 

leads to a order contribution, so it cannot enter here. Thus, we contract the (j) field from Eq. (j66p . using Eq. ([8]), 
with the LO mass coupling of (f) given by 

L(m)o = -mcj). (68) 
Again, we obtain a one-graviton exchange diagram, that which is depicted in Fig. 1(b). This diagram is evaluated by 

Fig. l{b) = [Si • vi X n + + [1 O 2]. (69) 



ov' v'Q- 



(a) 



(b) 



FIG. 1. LO spin-orbit interaction Feynman diagrams. The heavy solid lines represent the point particles worldlines. The oval 
gray and spherical black blobs represent the spin and mass couplings on the worldline, respectively. The solid and dashed lines 
represent the and Ai fields, respectively. All diagrams must be included with their mirror images. 
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Here again, we have the Fourier integral given by Eq. (jA2l) . One easily verifies this diagram also scales as v^. No 
other diagrams can be constructed at this order. Thus, the LO spin-orbit potential is of the 1.5PN order, and it 
follows from these two Feynman diagrams of one-graviton exchange. We point out that both diagrams here as well as 
all diagrams contributing to the spin-orbit interaction are always considered together with their mirror images. Thus, 
the [1 o 2] term will be omitted henceforth from all diagram values. 

The two diagrams here add up to the LO spin-orbit Lagrangian, which at this PN order is trivially Legendre 
transformed to a potential, namely VgQ = —L^q, with the momenta conjugate to the coordinate velocity being just 
the Newtonian momenta, i.e. p = mv. Thus, we obtain the well-known result for the LO spin-orbit potential 

^so = [Si • (vi X n - 2V2 X n) + S°'n'] + [1 O 2] = -^Si • Qvi x n - 2v2 x + [1 o 2], (70) 

where on the right-hand side the leading term of NW SSC from Eq. is substituted in, and the result is then 
in agreement with the canonical result found in e.g. Eq. (4.11a) of [32|. We would like to stress that unlike the 
LO spinl-spin2 potential, the LO spin-orbit potential already requires the use of SSC, a fact which demands our 
special attention. It is the 5°* spin entry that appears at the leading spin coupling which adds complication to our 
calculation. 

Let us consider then the LO spin-orbit potential given in Eq. ([70]) with the use of the covariant SSC of Eq. (|42p . 
instead of the NW SSC. Upon the use of the covariant SSC from Eq. (|47)) in Eq. ((TO)) , the LO spin-orbit potential 
takes the form 

Y^g = • (2vi X n - 2v2 X n) + [1 2] = -2^J^Si • v x n + [1 2], (71) 

where in the last equality we have used the relative velocity v = vi — V2. However, the known canonical result for 
the LO spin-orbit Lagrangian obtained with the use of covariant SSC takes the form [5l| 

L^g = 2^Si-vxn+iSi.vi xvi + [lo2], (72) 

where Va = dva/dt is the acceleration of the particle. It is easy to see that Eq. (|72|) differs from the result obtained 
via the EFT calculation in Eq. (j7ip by an additional acceleration dependent piece. One may be reluctant to have such 
terms, since they call for the use of the EOM at the level of the action, a procedure which is known to be incorrect 



in many (in)famous examples. However, as originally noted by 52|, and also treated by many others, e.g. |53l |54|. 
a substitution of low order EOM in higher order terms in the level of the Lagrangian is a correct procedure in GR, 
and is equivalent to making a coordinate transformation. Here, we may eliminate the acceleration terms using the 
LO EOM given by 

G(toi+TO2) 

V = a = ai - a2 = 5 n, (73) 



and we see that the canonical result obtained with the NW SSC in Eq. (|70|) is recovered. 

Alternatively then, one can recognize that the canonical result obtained with the covariant SSC in Eq. (|7^ . is 
related to the canonical result obtained with the NW SSC in Eq. (iTOt . by the following coordinate transformation at 
the level of the Lagrangian 

^ -S„ X Va ^ v„ ^ v„ + -^Sa X Va, (74) 



2ma 2rna 

where the spin precession S is considered at Newtonian order as 0. Equivalently, they are related by the following 
canonical transformation at the level of the Hamiltonian 

r- ^r- ^ -^Sa X Pa, (75) 



2m2 



Pa ^ Pa - 7r~Sa X Pa. (76) 
2mQ 

So now it becomes clear how to relate the EFT result obtained with the covariant SSC, to the canonical results 
obtained with the use of the covariant SSC or NW SSC. To relate Eq. (|7T|) with Eq. ([72]) one should make the 
following noncanonical change of variables at the level of the Hamiltonian: 

Pa ^ Pa ^" 7; X p^, (77) 

2ma 
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transforming only p, and not r, whereas to relate Eq. (j7ip with Eq. (I70p one should make the following noncanonical 
change of variables at the level of the Hamiltonian: 

fa ^ Ta + -i^-Sa X Pa, (78) 

transforming only r, and not p. The transformation in Eq. (|78|) is just similar to the flat spacetime mapping from the 
center of mass coordinate variable related with the covariant SSC, to that related with NW SSC 39]. To conclude, 
we have seen that already at the LO of the spin-orbit interaction the result obtained with an EFT calculation using 
the covariant SSC leads to a noncanonical result, related to the canonical results by noncanonical change of variables. 



V. NEXT-TO-LEADING ORDER SPIN-ORBIT INTERACTION 



In this section, we present the calculation of the Feynman diagrams, which contribute to the NLO spin-orbit 
interaction potential between the binary constituents. In the NLO spin-orbit potential, which is of the 2.5PN order, 
we have 0{G) and O(G^) contributions, so there are three diagram topologies contributing. First, we have the single 
topology at 0{G) of one-graviton exchange, and then we have the two topologies at O(G^) H^]. At 0{G'^) we first 
consider the "V" topology, i.e. diagrams of two-graviton exchange. Then, we must consider the "Y" topology, namely, 
diagrams that contain a three-graviton vertex. Note that after stripping off worldlines from our diagrams, we are 
always left with connected diagrams, and no graviton loops. We use the Feynman rules presented in Sees. HIl and IIIII 
for the propagators, propagator correction vertices, three-graviton vertices, and the worldline mass and spin couplings. 
In addition, we have to keep track in each diagram of the possible numerous contractions, and the symmetry factors 
arising from the worldline couplings at 0{G^) ^55]]. We evaluate each diagram in turn, starting from the simplest 
one-graviton exchange diagrams and ending with the more complex time dependent cubic self-gravitational interaction 
diagrams. 



A. One-graviton exchange 

For the NLO of the spin-orbit interaction, we start by still considering one-graviton exchange diagrams. Calculating 
now to order (2.5PN), there are seven diagrams to be evaluated as shown in Fig. 2. 

We should include diagrams with a spin coupling of -/j(s)2 from Eq. and the subleading mass coupling from 

Eq. (gll) given by 

im(3) = ]^mA,v'v^. (79) 

We obtain the diagram shown in Fig. 2(al). The value of this diagram is given by 

i^zg. 2(al) = -^^Si-V2Xnw2. (80) 

The only integral in k which arises is a Fourier transform, evaluated using the formula in Eq. (|A2p . This is the case 
for all one-graviton exchange diagrams that do not include propagator corrections, i.e. for Figs. 2(al)-2(a4). 

For the NLO of the SO interaction, we should also include diagrams with a spin coupling of L(s)3 in Eq. ([66)1 and 
the mass couplings of (j) given by 

3 

L(ra)2 ^ -l^mcjiv'^ , (81) 

and of aij from Eq. (|67| . Contracting each of the fields </> and aij, we obtain the diagrams shown in Figs. 2(a2.1) and 
2(a2.2), respectively. The values of these diagrams are given by 

Fig. 2(a2.1) = ^^^^ [Si • vj x n + S^W] wf, (82) 
IGth 

Fig. 2(a2.2) = ^Si • [vi x n w| - V2 x n (vi • V2)] . (83) 



15 



Next, we must consider the spin coupling from Eq. (j57p at order given by 

i(s)4 - -\s'^'d,A,v' + i^^aovl,. (84) 

Note that both terms here contain temporal spin entries. Moreover, note the second term here that first contains time 
dependence. We should then also include diagrams with the spin coupling of i(s)4 together with the mass coupling 
L[m)i of Eq. (|64)) . Such a diagram is depicted in Fig. 2(a3) and is evaluated by 

Fzg. 2(a3) = [^?V^(vi • v^) + Sfv^y, ■ n)] + '^-^dtSfv^. (85) 



Here, when contracting the d^Ai term of Eq. (|84p. we make an integration by parts to drop the time derivative from 
the delta function on the Fourier transform integral, and the temporal spin component Soi. Only then do we perform 
the integration over time, and use the delta function. Note that the time derivative of the temporal spin component 
may yield an acceleration dependent term from the implicit velocity power it carries, see e.g. the SSC in Eq. (|47p . as 
well as a spin precession term. Here, we may eliminate the acceleration terms using the LO EOM given in Eq. (j73p . 
which will result in a nonlinear O(G^) contribution. By such a substitution we will implicitly have a coordinate 
transformation made. As noted, the time derivative may also yield a spin precession term that requires the use of the 
LO spin EOM. However, the LO spin EOMs scale as [3^ 

dS Gm 



dt 



Sv. (86) 



This will yield a subleading order contribution and should therefore be dropped here [ISj . Because of the term with 
the time derivative this diagram may be evaluated in two ways. The contraction with the doAi term may be evaluated 
alternatively by flipping the time derivative between the two particles, namely, by using the identity 

dhdt2 dt,S{h - t2)f{ti)g{t2) = - j dtidt2 dtJih - t2)f{h)g{t2). (87) 



u a 



(al) 



(a2.1) (a2.2) 



(a3) 



(a4) 



(bi) 



(b2) 



FIG. 2. NLO spin-orbit interaction Feynman diagrams of one-graviton exchange. The double line represents the (Jij field. The 
encircled cross vertex corresponds to a propagator correction. All diagrams should be included with their mirror images. 
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Evaluating the diagram this way yields the following seemingly different value for the diagram in Fig. 2(a3), given by 

Fzg. 2{aiUt. = [STn^i^i ■ ^2) + S"Mi^2 ■ n)] - ^^S'M- (88) 

Hence, we conclude that an acceleration dependent term and the use of EOM are inevitably involved in evaluating 
this diagram (though we may avoid a spin precession term on the later evaluation). We stress that the two evaluations 
can be related by a total time derivative, and are thus physically equivalent. 

Next, considering spin couplings at order v^, we take into account the flat spacetime subleading substitution in the 
SSC of Eq. (|48p . It may make the S*"* term of Eq. ((58)) become a spin coupling of order given by 

i(s)5 = (89) 

Contracting this spin coupling, using Eq. PH)) . with the LO mass coupling of Eq. we may get the diagram shown 
in Fig. 2(a4), which is evaluated by 

Fig. 2(a4) = ^^S°'n\ (90) 

We should stress that we do not have this contribution, if we use the covariant SSC at the level of the action. 

There are further contributions to the NLO SO interaction from the one-graviton exchange sector, which arise from 
corrections to the potential graviton propagators. These are the self-gravitational vertices of Eq. ([T^ . containing 
two time derivatives, and hence suppressed by v'^. First, we consider diagrams that contain the LO spin coupling of 
Eq. ([63]) . contracted with the propagator correction in Eq. (fTSj) . and the mass coupling of Eq. ([64|. using Eq. ([9]). 
Such a diagram, which is shown in Fig. 2(bl), is evaluated by 

Fig. 2(61) = ^j^Si • [-V2 X n (vi • V2) + vi x V2 (v2 • n) + 3v2 x n (vi • n)(v2 • n)] + ^^Si • vi x a2. (91) 

Note that here we have to consider two similar contractions. To compute this diagram we first have to use the identity 
of Eq. (|87p twice - for each time derivative, so that we get dt^ and dt2 ■ Then, we make two integrations by parts to 
drop the time derivatives on the respective Fourier integrals, or on or V2, respectively. Then, we apply the time 
derivatives. Thus, we get a spin precession term that drops on the use of the LO spin EOM from Eq. ([86| . and an 
acceleration dependent term with a2, to be substituted using Eq. (f73| . as explained above. Only after using the time 
derivatives do we perform the time integrations, using the delta functions. Finally, we arrive at a Fourier integral, 
whose value is given in Eq. (|A3[) . 

As for the diagrams including the propagator correction of Eq. (|14p , they contain the spin coupling of 4> in Eq. (|66|) , 
and the mass coupling of Eq. all contracted using Eq. Such a diagram, shown in Fig. 2(b2), is evaluated by 



— ^Si • vi x n (vi • V2) + vi X V2 (vi • n) - 3vi X n (vi • n)(v2 • n) + — — ! 
2r^ 2r 



[Sfn\w, ■ V2) + 5?'«^(vi • n) - 35«n'(vi • n)(v2 • n)] - 9^dtSfv\. 



Here also we have two similar contractions. The evaluation is similar to that of the previous diagram, except that 
the time derivative falls on a temporal spin component S^"^, which may yield an additional acceleration term upon 
the use of SSC, e.g. in Eq. (|47)) . as in the evaluation of Fig. 2(a3). In addition, we get terms containing S'l^v\ and 
dtSi^n^, which can be seen from power counting to eventually drop at the implementation of the SSC. Here again, 
acceleration dependent terms and the use of EOM are inevitable. 



B. Two-graviton exchange 

At NLO in the SO interaction nonlinear contributions must also be included. At 0{G^), we consider first the 
two-graviton exchange, i.e. the diagrams of the V topology. There are four such diagrams to be evaluated here as 
shown in Fig. 3. They arise from two-graviton couplings of both mass and spin. First, we consider the two-graviton 
mass couplings, which are simpler. At order we have the two-graviton mass coupling from Eq. ([27]), i.e. 



L{m)3 = ■m(t>Aiv\ 



(93) 
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(al) (a2) 




FIG. 3. Nonlinear NLO spin-orbit interaction Feynman diagrams of two-graviton exchange. These diagrams should be included 
together with their mirror images. 



This must be contracted with the LO spin coupling of Eq. (IMl) . using Eq. ©, and the LO mass coupling of Eq. 
using Eq. ([8]), to obtain the diagram depicted in Fig. 3(al), which equals 

Fz5.3(al) = ^^^^Si.V2Xn. (94) 

The calculation of the two-graviton exchange diagrams is rather simple since the two momentum integrations just 
factorize into the two simple Fourier transforms of the one-graviton exchange, i.e. that of Eq. (IA2p . 
We must also take into account the two-graviton mass coupling of Eq. (|26)) given by 



i(m)2 = -Y^^- (95) 

This should be contracted with the spin coupling of (j) in Eq. (j66| , and LO mass coupling of Eq. ([68| , using Eq. ([8|) , 
to obtain the diagram depicted in Fig. 3(a2), which equals 

Fig. 3(a2) = -^^^^^^ [Si • vi x n + S°,'n'] . (96) 

Note that here we have two similar contractions. 

Next, we should consider the two-graviton spin couplings. First, we have from Eq. (|60p the LO two-graviton spin 
coupling given by 

L^s)4 = S''F,jcb-^S'^A,djcl). (97) 

The spin coupling of Eq. ()97p scales as w^, so diagrams containing it must be contracted with the mass couplings of 
Eqs. (j68| and (|64l). using Eqs. (jS]) and ([9]). Such a diagram is shown in Fig. 3(bl), and it equals 

lOG^TO^ 

Fig. 3(&1) = ^3-^ Si • V2 X n. (98) 

In addition, we have a contribution from the 0{v^) two-graviton spin coupling of Eq. (j6ip . given by 

i(s)5 = 25°'9,#. (99) 

Note that this two-graviton spin coupling appears with the S*"* entry, so the use of SSC is involved. This two-graviton 
spin coupling should be contracted with two of the mass couplings of Eq. (|68p , using Eq. ([5]) . This forms the diagram 
depicted in Fig. 3(b2) and given by the value 

Ftg. 3(62) = ^^1^ Sfn\ (100) 
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Here, the symmetry factor is ^ but we have two similar contractions. 

Recall that there are metric field corrections in the temporal 5*'° entries from Eq. (l62t . These may be considered 
as additional contributions to the two-graviton spin couplings in Eqs. ([97|) and (f99|) here, arising from the S'^'^dicj) 
coupling in Eq. ([55| . if the SSC is applied at the level of the action. There is an additional contribution to the SAcj) 
couphng in Eq. ((97)) . which was already found in the NLO spinl-spin2 sector [l^, and to the Scp"^ coupling in Eq. ([99)) . 
which did not appear in the NLO spinl-spin2 sector. We will consider these additional contributions in Sec. IVH 



C. Cubic self-gravitational interaction 

Finally, we must include the nonlinear contribution arising from diagrams, where a three-graviton vertex mediates 
the interaction, i.e. the second topology at O(G^) - the Y topology [IS]- As shown in Fig. 4, there are seven 
such diagrams to evaluate in the NLO SO interaction. All of these diagrams contain one-loop integrals and require 
regularization. The regularization can be understood to arise from the fact that unlike n-graviton exchange between 
two worldlincs, in which the graviton momenta have a maximal cutoff, set by the reciprocal of the separation distance 
r between the two worldlines, here it is possible for the three-graviton vertex to be arbitrarily close to the worldlines 
(in position space, of course), thereby yielding arbitrarily large momenta for the gravitons in the cubic interaction. 
It may also be interpreted so that these diagrams have a dressed vertex, namely, that the two vertices of the same 
object can be replaced by a single dressed vertex [56|. The loop integrals are handled with dimensional regularization 
by the usual techniques, see e.g. 113, and Appendix El for the required formulas. It should be stressed that the NRG 
field parametrization we employ here certainly makes the computation of the cubic gravitational interaction diagrams 
feasible by hand, compared with the long calculation required in e.g. [l2| or [l3|, which relies on the aid of automated 
computations. 

The LO three-graviton vertices of Eq. (fTB)) scale as [l3|. First, we consider the (pF"^ vertex of Eq. (US]). This 
vertex should be contracted with the LO spin coupling of Eq. and the mass couplings of Eqs. and (|64l) . 

using Eqs. ([8]) and ([9]) respectively. However, one should consider all possible insertions of the mass couplings on 
the two worldlines. One can have both mass couplings on the same worldline as shown in Fig. 4(al), or each mass 




(al) (a2) (a3) 




(bl) (b2) 




(cl) (c2) 



FIG. 4. Nonlinear NLO spin-orbit interaction Feynman diagrams with a three-graviton vertex. The gray encircled cross vertex 
corresponds to the time dependent cubic vertex. These diagrams should be included together with their mirror images. 
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coupling on another worldline, where in the later possibility an additional permutation can be made, thus yielding 
the two diagrams appearing in Figs. 4(a2) and 4(a3). These diagrams are evaluated, respectively, as 

2 

Fig. 4(al) = -^^^ Si • V2 x n, (101) 
F^g. 4(a2) = -^^^^ Si • x n, (102) 

Fig. 4(a3) = Si • vi X n. (103) 

In each of these diagrams, we have two similar contractions. In all of them we have to evaluate a momentum integral 
over two momenta. To calculate it, we first perform the integration that corresponds to a one-loop integral and then 
handle the remaining Fourier transform integral over the other momentum. In Fig. 4(al), we use the one-loop integral 
given in Eq. (|A6[) . and then we should use the Fourier integral of Eq. (|A4p . We note also that in Fig. 4(al) the 
contraction of the second term in the vertex of Eq. (jl8p cancels out, due to the antisymmetry of the spin tensor. In 
Figs. 4(a2) and 4(a3), we use the one-loop integral given in Eq. (jA7|) and again the Fourier integral of Eq. (IA4|) . 

Next, we also have the cr0^ vertex given in Eq. (1191) . This vertex is contracted with the u term of the spin coupling 
in Eq. (1551) . using Eq. ([TU|) . and with two of the mass couplings in Eq. using Eq. Here also, one can either 
have both mass couplings on the same worldline, leading to the diagram shown in Fig. 4(bl), or each mass coupling 
on another worldline, leading to the diagram shown in Fig. 4(b2). These diagrams are evaluated, respectively, as 

i^*5-4(61) = -^^ Si-vi xn, (104) 

Fz5.4(62) = -^l^^^Si.vixn. (105) 

The symmetry factor of Fig. 4(bl) is ^. Both diagrams have two similar contractions. In both diagrams one is 
required to use both one- loop integrals given in Eqs. (IA6|) and (|A7[) . and the Fourier integrals given in Eqs. (jA4[) and 

(US). 

At the NLO SO interaction we have in addition a contribution from another three-graviton vertex that contains a 
single time derivative. This is the Acj)^ vertex given in Eq. (|20|) . which scales as v'^ . This vertex is contracted with the 
LO spin coupling of Eq. using Eq. ([3]), and with two of the LO mass couplings of Eq. using Eq. ([S]). Here 
again, just like the diagrams that include the a4>^ vertex shown in Fig. 4(b), one can either have both mass couplings 
on the same worldline, leading to the diagram shown in Fig. 4(cl), or each mass coupling on another worldline, leading 
to the diagram shown in Fig. 4(c2). These diagrams are evaluated, respectively, as 

Fig. 4(cl) = Si • V2 X n, (106) 

Fig. 4(c2) = Si . VI X „. (107) 

Here, just like with the propagator corrections in Figs. 2(bl) and 2(b2), we have the further complication of the time 
derivative on the gravitational vertex. Just like we did there, we first use the identity of Eq. (j87p and then integrate by 
parts to drop the time derivatives on the Fourier integrals. Subsequently, we apply the time derivatives and only then 
do we perform the time integrations, using the delta functions. In both diagrams one uses both one-loop integrals 
given in Eqs. (jA6l) and (jA7p . and the Fourier integrals given in Eqs. (jA4|) and (jA5[) . Figure 4(cl) has a symmetry 
factor of i and two similar contractions. Figure 4(c2) has two different contractions, where the second contraction, 
in which one contracts the time derivative with the mass coupling of particle 2, vanishes, due to the antisymmetry of 
the spin tensor. 

We note that in the NRG field parametrization used here the cubic gravitational interaction sector does not contain 
temporal spin entries, which is advantageous in this more computationally complex sector. Finally, it is also worth 
noting that had we used the fully harmonic gauge rather than the linearized harmonic gauge, which would affect the 
three-graviton vertices used here, we would still obtain the same values for the cubic self-gravitational interaction 
diagrams. Hence, we can conclude that the total NLO spin-orbit interaction result that we get would be unchanged. 



VI. NEXT-TO-LEADING ORDER SPIN-ORBIT POTENTIAL 



Summing up all of the contributions from the Feynman diagrams, we obtain the NLO spin-orbit interaction La- 
grangian. As was noted already, we can consider our Lagrangian as a Routhian, since the Legendre transformation 
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to a Routhian with respect to the spin conjugate coordinate is trivial. Here, we present only the Si dependent part, 
and one should add a similar term for S2 by exchanging particle labels 1 and 2 in the following expression. Our NLO 
spin-orbit interaction Lagrangian for a binary system is then given by 
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^so 



Gm2 



Gm2 



Sr 



Sr 



vi X n 



1 



1 



-Vi • V2 



2 
V2 



-(vi • n)(v2 • n) 



V2 X n Vi • V2 



3(vi • n)(v2 • n; 



vi X V2 1 -vi • n + V2 • n 
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-V2 X ai 



^G2m2_ 
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vi X n ^ 
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V2 X n 



:m2 
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G^rn2 „o 



STn' [mi + 2m2] , 



(108) 



Here, we give the result with the acceleration dependent terms, and we do not substitute the SSC in the 5"°' dependent 
terms, so that the last two lines of Eq. (llOSp represent the SSC dependent part. Naturally, the additional contributions 
of the field corrections in the 5'°' entries from Eq. ([5^. which may be considered in the two-graviton exchange sector, 
if the SSC is considered on the level of the action, are also not included in this form. 

First, we substitute in Eq. (|108p the accelerations using the LO EOM given in Eq. (|73p as was discussed already. 
We obtain 
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2 2 
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(109) 



This result is to be compared with Eq. (4.11b) of [33 |. 

The NLO SO Lagrangian should be Legendre transformed with respect to the velocities Vq in order to obtain a NLO 
SO Hamiltonian that we can compare. We recall that the Legendre transformation with respect to the coordinate 
conjugate to spin is trivial and may be considered as done once the spin is treated as an independent variable and 
the Lagrangian is taken with the opposite sign. The part of the Lagrangian required for the obtainment of the NLO 
SO Hamihonian iJ^o ° reads 



L — Ln 

where Ln is the Newtonian Lagrangian given by 



T ^rLO , tNLO 

^EIH — VsO + ^'^ 



^SO ' 



(110) 



1 ^ 



Gmim2 



(111) 



and Leih is the Einstein- Infeld- Hoffmann Lagrangian, which is the IPN correction of the Lagrangian [l2|, given by 



L 



EIH — 7: 



a=l 



GmiTO2 
2^ 



[3vl + 3vl - 7vi • V2 - (vi • n)(v2 • n)] - 



G^mim2(mi -I- m2) 
2^2 ■' 



(112) 



and where V^q and L^q'-' are the LO and the NLO spin-orbit couplings, which are given here in Eqs. (|70p and (|109p 
respectively. The canonical momenta conjugate to the coordinate velocities are given by 



Pa 



dL 

dWa ' 



(113) 
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so that we have 



pi = mivi 



-Vl 



Gmim2 



3vi - -V2 



V2 ■ n \ G 

n (m2Si X n + 2miS2 x n) + 



arNLO 

9vi 



(114) 



where we note that here 5**° is still considered as an independent degree of freedom, and P2 is obtained from pi by 
exchanging particle labels 1 -o- 2 in the above expression. By inverting the relation between conjugate momentum to 
the coordinate velocity to the desired PN order, we have for the coordinate velocity the following PN expansion in 
terms of the canonical momenta 1321 
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(115) 



This can be obtained equivalently by using vi = with the relevant Hamiltonian. Again, V2 is obtained from Vi by 
exchanging particle labels 1 2 in the above expression. Finally, we perform the Legendre transformation, namely, 
we substitute in the coordinate velocities using Eq. (jlisp in the following expression: 



H 



^ Va • Pa - ^■ 



(116) 



a=l,2 



Note that the a?P terms get canceled and do not contribute here. 



dVa , , 

The obtained Hamiltonian then takes the following form: 

H = Ho + Hso 



H: 



ss, 



(117) 



where Ho is part of the orbital interaction Hamiltonian, Hso is part of the spin-orbit (linear in spins) Hamiltonian, 
and Hss is part of the spin-spin (quadratic in spins) sector of the Hamiltonian. Thus, the NLO spin-orbit Hamiltonian 
we obtain equals 
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Note the first term appearing in the SSC dependent part of the Hamiltonian, arising from the LO SO interaction. We 
need to pay special attention and keep track of this 5*°* term, as it contributes to the NLO SO interaction too. 

Now, we are going to eliminate the S^^ entries from the Hamiltonian using the covariant SSC in Eq. (j62p . utilizing 
Eq. (jll5|) . and the metric for spinning binary black holes in harmonic coordinates, e.g. as in (58| . It is well known 
that this will yield noncanonical variables in the Hamiltonian. However, we can employ noncanonical transformations 
of variables that relate the variables to canonical ones [l^. Therefore, after the elimination with the covariant SSC, 
we obtain the following Hamiltonian with noncanonical variables: 
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We go on to apply the noncanonical transformations of variables that read 
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(120) 
(121) 



These transformations are just the generalization of the mapping between covariant and NW SSC variables in flat 
spacetime [s^. The transformation of the spin variable in Eq. p20p is just similar to that of flat spacetime, whereas 
for the center of mass coordinate in Eq. (|121l) a higher order PN transformation, than what we had in Eq. (|78l) . 
generalized for curved spacetime, is required [23, 32]. We argue that now we have arrived at a canonical Hamiltonian. 
As was mentioned before, we compare our result with that given by Eq. (4.11b) of ^32*], or in our notations 
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The difference between the part linear in Si in H^q^ and that of the spin-orbit Hamiltonian obtained here in Eq. (|119p . 
after the application of the noncanonical transformations in Eqs. (|120p and (|12ip . i.e. 
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is given by 
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If the Hamiltonian obtained here is physically equivalent to that of [32], i.e. if it yields the same EOM, there should 
exist an infinitesimal generator g of a canonical transformation (CT) such that 



dg 
dt' 



(125) 



So let us construct a suitable infinitesimal generator following the discussion of PN canonical transformations in 
Appendix |B] . Such a generator must be a scalar, linear in the spin vector, and constructed out of the spin vector, 
Pi, P2, and n. Since the spin already carries in it, we note that there may be three types of generators: ©(G^p^), 
O(G^p^), and 0{G^p^). However, only the 0{G^p^) type contributes appropriate terms. Hence, we consider the 
following ansatz for the O(G^p^) generator: 



G7772 _ 
g = Si 



31- 



Pi X p2 , Pi X n 



7)7 1 7772 



7771 



52- 



Pi n 

7771 



.93- 



P2 ■ n 

7772 



P2 X n 

7772 



Pi • n p2 • n 

54 1- 95 



mi 



7)72 



(126) 



This generates both 0(G) and O(G^) terms of the form that appears in the NLO spin-orbit sector: 0(G) terms from 
differentiating the coordinates with respect to time, and O(G^) terms from differentiating the momenta with respect 
to time, or more precisely p terms, which become O(G^) terms upon the use of LO EOM. For example, from the gi 
part of the generator in Eq. (|126p . we get the following transformations for the canonical variables ri and pi: 



ri 



ri -51- 



G7772 Si X P2 



7771 7)72 
G7772 Pi X P2 

Pi Pi + 91 — 5- Si • 

777 1 7772 



(127) 



23 



Thus, we plug in Eq. (jl25l) our ansatz for g from Eq. (I126p . and we compare that to Eq. (|124[) . Comparing 0{G) 
terms gives 

51 = ^, 32 = 0, ff3 = -^, 54 = 0, 55 = 0. (128) 

This ehminates all of the 0{G) terms, as well as the 0{G^) term from the difference in Eq. (|124p . i.e. finally 

AH^^" ^ (129) 

in agreement with the canonical result in (32j . 

It is interesting to note that had we used the covariant SSC at the level of the action, rather than considering 
the temporal spin entries S''^ as independent degrees of freedom until their elimination at the Hamiltonian level, we 
would have arrived at the same Hamiltonian appearing in Eq. (|119p . We also note for completeness that the physical 
equivalence of our result with that of (32j or 33 1 can also be shown by computing and recovering the EOM, using Dirac 
brackets. However, working with variable transformations is more efficient, which becomes crucial when computing 
similar or higher order spin corrections. 



VII. CONCLUSIONS 



In this paper we applied an EFT approach to calculate the NLO gravitational spin-orbit interaction between two 
spinning compact objects. The NLO spin-orbit interaction was the last NLO conservative spin correction to be 
computed within the EFT approach, due to its higher complexity with respect to the NLO spinl-spin2 and spin- 
squared sectors, previously computed for the first time within the EFT approach. The calculation of the NLO 
SO sector involves the evaluation of 18 Feynman diagrams, 7 of which correspond to one-loop diagrams, including 
nonstationary cubic self-gravitational interaction, as well as other nonstationary contributions. It is also worth noting 
that in the NLO SO sector one cannot avoid having acceleration dependent terms, and using EOM to eliminate 
them. The NRG field decomposition, first applied in the computation of the NLO spinl-spin2 interaction, in terms of 
which the EFT calculation is carried out here, facilitates the calculation considerably. We recall that the NRG field 
decomposition in the harmonic gauge has many advantages, such as simple propagators for the scalar and vector fields, 
which dominate in the SO interaction, and the vanishing of mixed 2-point functions. In particular, the NRG field 
decomposition makes the treatment of the cubic self-gravitational interaction relatively simple, so that there is no need 
to rely on automated computations at all. In addition, the NRG field decomposition illustrates the coupling hierarchy 
of the different gravitational field components to the spin and mass sources, which is helpful in the construction of 
Feynman diagrams. 

However, the main obstacle in the calculation of the PN spin corrections is the treatment of SSC. Therefore, the 
fact that, unlike the LO spinl-spin2 or spin-squared effects the LO spin-orbit effect already involves the use of SSC 
presents a major complication in this sector. We showed here that already at the LO SO sector, this results in 
ambiguities related with different choices of SSC, which can be resolved in terms of noncanonical changes of variables. 
Moreover, the NLO SO sector requires the application of SSC at higher orders, whereas SSCs at higher orders were 
not considered until recently, and they should be further explored. Finally, we have shown here explicitly how to 
relate the EFT derived spin results to the canonical results obtained with the ADM Hamiltonian formalism. This is 
done, similarly to the LO SO sector, by using noncanonical transformations, as well as canonical transformations at 
the level of the Hamiltonian. We note, however, that it is more difficult to apply on the Hamiltonian noncanonical 
transformations compared to canonical ones. Yet, working with transformations of variables is more efficient than 
resorting to the EOM and Dirac brackets, and that becomes crucial when unknown higher order spin corrections are 
to be approached in the future. 
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Appendix A: Dimensional regularization and loop integrals 

Throughout the computation of the contributing Feynman diagrams, we encounter two types of momentum integrals 
that need to be evaluated: Fourier integrals that arise from the Fourier transforms of the propagators, and one-loop 
integrals, which arise from the cubic gravitational interaction. Both types of integrals are evaluated using dimensional 
regularization [H^l- In order to evaluate the Fourier integrals, one should use the d-dimensional master formula given 
by 
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This formula can easily be derived using Schwinger parameters 57] . From this master formula, we obtain the following 
required Fourier integrals: 
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The loop integrals required for the computation of the Feynman diagrams that contain thrce-graviton vertices are 
given by 

d^k k^ ^^^^ 
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These integrals can be obtained using both Feynman and Schwinger parameters [5 
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Appendix B: PN canonical transformations 

Let us elaborate on the statement made by Eq. (|125p and derive it here. Recall that our Lagrangian is acceleration 
dependent, due to its high PN order (2PN order and beyond). We know that for a Lagrangian of the form L = L{q, q, t), 
i.e. a standard Lagrangian, which depends on some coordinate q and velocity q, a total time derivative of an arbitrary 
function F{q,t), which depends on the coordinate q and time only, can be added to the Lagrangian L, such that 
L{q,q,t) + ■^F{q,t) is physically equivalent to L, i.e. it yields the same EOM as L. Similarly, it is easy to show that 
for an acceleration dependent Lagrangian of the form L = L(q, q, q,t), a total time derivative of an arbitrary function 
F{q, q, t), which depends also on the velocity g, can be added to the Lagrangian L, such that L{q, q, q, t) + ^F{q, q, t) 
yields the same EOM as L. 

An addition of such a total time derivative to the Lagrangian corresponds to applying the infinitesimal CTs given 

by 

q^ q + da^g{q,p), 
d 

p^p-da—g{q,p), (Bl) 

on the canonical variables (9,^) of the Hamiltonian, where da is an infinitesimal of the transformation parameter, and 
g — g{q,p) is the infinitesimal generator of the transformations. If one would like to make an addition to the ^PN 
part of the Hamiltonian, g should be of the same PN order, or more precisely it should scale as v^t. The generator 
g should be given in terms of our canonical variables, and it should not depend on time explicitly, since our PN 
Hamiltonian does not. It can be directly checked that the mappings in Eq. (jBll) describe a canonical transformation 
from the so-called "direct conditions" for restricted canonical transformations [46J, i.e. CTs that do not contain time 
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explicitly. Note that since g '-^ , the canonical momenta p that will appear in g will effectively be the Newtonian 
momenta, i.e. p = mv. 

Now, let A be some function of the canonical variables, which does not depend on time explicitly. Then, it follows 
that dA — da{A, g} [i^, so that the change in A due to a finite CT is given by 

2 3 

AA = a{A, g} + ^{{A, g}, g} + ^{{{A, g}, ff}, ff} + • • • . (B2) 

Since g scales as u", all orders beyond linear in a do not contribute to the considered PN order, so to that order we 
have AA = a{A,g}. From this we find that the change in the Hamiltonian due to the finite CTs, which are of the 
form 

dg 

q+a—, 
op 



is given by 



p^p-a^, (B3) 



AH = a{H,g}^~a^, (B4) 



since our H and g does not depend on time explicitly. Thus, we have arrived at Eq. (|125p . up to a numerical coefficient, 
which is the finite transformation parameter. 
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